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Abstract 

Understanding the charge carrier transport in the disordered organic molecular semiconductors 
is a fascinating and still unresolved problem in modern condensed-matter physics, yet has an 
extremely important bearing on their application in "plastic" based field effect transistor, light 
emitting diodes and lasers. The most contentious issue in this subject is the nature of the charge 
carriers i.e. whether they are bare electrons or dressed electrons, known as polarons. Here we 
show for the first time, by means of simple experiments that polaron moving in static disorder is 
responsible for charge carrier transport in disordered molecular semiconductors. 

PACS numbers: 72.80.Le, 72.15.Cz, 72.20.Ee 
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Since the discovery of "plastic" based transistor [l|, Q and light emitting devices 0, 
disordered molecular semiconductors are being widely investigated. These materials are 
composed of different types organic molecules starting from very small molecules to very 
large biomolecules like proteins and DNA, held together loosely by weak van der Waals 
type intermolecular force. The absence of long range order in these materials leads to 
the localization of the electronic wave function and results transport of charge carrier via 
thermally activated hops with carrier mobility fi as a strong function of temperature and 
electric field F following a universally 6] observed Poole- Frenkel(PF) behavior fi(F, T) = 
/i(0, T) exp(7\/F), where /i(0,T) is the temperature dependent zero field mobility and 7 
is field activation of the mobility and recently, in light of theoretical studies jl, 0, E| we 
have shown [lo| another important feature i.e. (iv) charge carriers are spatially correlated in 
these materials. Considerable debate is going on the fundamental question like, what is the 
nature of charge carriers in disordered molecular semiconductors'! Answer to this question 
will be a key input for understanding the charge transport mechanism and the development 
of new materials and new devices based on these materials. An indication of the difficulty 
in finding the nature of charge carriers and understanding the transport processes in these 
materials is evidenced by the conflicting viewpoints. At present, there are three viewpoints 
on the nature of charge carriers, which are (i) dressed electrons (known as polarons) due to 
strong electron-phonon interaction and the temperature dependence of the /1 arises due to 
polaron binding energy (if) bare electrons and temperature dependence of \i is due to 
energetic and spatial disorder [lj] and (iii) third viewpoint 0, [l^ is the combination of these 
two opposite viewpoints and the charge carriers are polarons moving in disorder and the 
temperature dependence of \i is due to both static disorder and polaron binding energy. It has 
been shown that polaron based models results unacceptable values of polaron binding 
energy and strong nearest neighbor intermolecular coupling if disorder is not included. As 
mentioned earlier, the intermolecular coupling in these disordered molecular materials is 
very weak and can't account for strong nearest neighbor coupling. If it is assumed that 
energetic and spatial disorder are responsible for thermal activation of \x and bare electrons, 
not polarons, are responsible for charge carrier transport, the probability of hopping is 
drastically reduced and again unreasonably strong intermolecular coupling is required to 
explain the experimental data. Here, we show that there is a natural solution to this apparent 
paradox and polaron is the nature of charge carriers, emphasizing the important role of both 
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disorder and electron-phonon interaction on the charge carrier transport in these disordered 
materials. 

We have chosen one of the most important (3, Q| molecular solid based on Copper ph- 
thalocyanine(CuPc) molecules because of their high chemical stability and ease to prepare 
the devices. Details about the CuPc-based single layer devices are given in Ref.Q]. Here, 
we report the experimental investigation on charge carrier transport in hole only devices 
based on met al/CuPc /metal structures. In this case, by properly choosing contacting met- 
als, current injection and transport due to holes have been investigated. Fig.l shows the 
temperature dependent current-voltage(J-V) characteristics in ITO/CuPc/Al structure. At 
low bias, J-V characteristic follow the slope of two in log(J)-log(V) plot, but as the bias 
increases the slope increases gradually to higher value due to space charge limited transport 
process|3], because the ionization potential of MePc(4.8eV) [3] is close to the work func- 
tion of ITO(4.75eV). A numerical workout 16] has been used by solving Poission's equation 
dF(x)/dx = ep(x)/e, current density relation J(x) = ep(x)n(x)F(x) and the PF relation 
of carrier mobility fi(F,T) = //(0, T) exp(7vF), simultaneously. Here, e is the dielectric 
constant of CuPc and p(x) is the hole density at position x and the boundary condition 
is taken }l5. at the ITO/CuPc interface. The simulated results are shown as solid lines 
in Fig.l at different temperatures. Inset of Fig.l shows the temperature dependence of 



conductivity a 
range hopping 



Tj for CuPc. Linear dependence of ln(a) on 1/T 1 / 4 validates the variable 
171 ]. cr(T) = a exp{ — (Tq/T) 1 / 4 }, where a is the temperature independent 
conductivity and To is the characteristic temperature. Hence the excellent agreement of the 
simulated results with experimental data presented in Fig.l establish the carrier transport 
process, which is thermally activated hopping with field dependent mobility. 

There are basically three models to explain the charge carrier transport in these ma- 
terials. In a very first approach, Gill 18] attempted to describe the universally observed 
thermally activated and PF behavior of \i using a phenomenological non-Gaussian disorder 
model(NGDM), fi(T,F — > 0) oc exp(— A/ksT), where A is the temperature independent 
activation energy. This model has been criticized [3, Q Q] for lacking a first principle 
explanation of PF behavior of fi and non-inclusion of another universal feature i.e. the 
Gaussian distribution of energy levels, g(E) oc exp (— E 2 /4a 2 ), where o is the r.m.s de- 
viation of hopping site energies. Bassler proposed 19| the uncorrelated Gaussian disorder 
model(UGDM), which describes the carrier transport as a biased random walk among the 
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dopant molecules with Gaussian-distributed random site energies due to positional and en- 
ergetic disorder. Although UGDM explains some features of experimental data and provides 
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support for PF behavior of carrier mobility, several discrepancies 
lated description of Gaussian disorder model. Garstein and Conwell 
spatially correlated potential is required for the description of PF behavior of mobility and 
subsequently Novikov et al.|7j have proposed a correlated Gaussian disorder model(CGDM) 
to describe the charge carrier transport in these materials. Recently, we have shown 10(] 
that the CGDM successfully explains the experimental data emphasizing, whatever be the 
nature of charge carriers, charge transport occurs by correlated hopping among the molec- 
ular sites distributed in Gaussian density of states(GDOS) of highest occupied molecular 
orbital(HOMO) and/or lowest unoccupied molecular orbital(LUMO) in these disordered ma- 
terials. In case of both Gaussian disorder models, the temperature dependence of zero field 
mobility is given by /i(T, F — > 0) oc exp (-Aa 2 /k 2 B T 2 ), where A is 2/3 in UGDmQ and 
3/5 in CGDmQ. The Aerrhenious plots result an apparent activation energy E c , which is 
E c = -k B {d(lnfi)/d(l/T)} = A in NGDm[l£| and Ba 2 /ksT, in case of Gaussian disorder 
models and B is 8/9^ in UGDmQ and 18/25^ in CGDmQ. If we take the contri- 
butions both from the polaron binding energy E po i, and from the static disorder E^, E c 
should be E c = E po i + Edis- Fig. 2 shows the temperature dependence of zero field mobility 
and it is clear that both ln(n) vs. 1/T and ln(fi) vs. 1/T 2 result nonlinear Arrehenious 
plots and it is impossible to distinguish between E po i and E^s and their respective values 
from E c . Here, we show that the interplay between these two contributions on the car- 
rier transport can be explored by studying the temperature dependence of /i over extended 
range of temperature and measuring atleast one contribution in E c , independently. It is 
clear from Fig. 2 that both, NGDM and UGDM or CGDM agree well with experimental 
data in higher temperature region(T > 140i^) and gradual change in slope in the Arrehe- 
nious plots can be explained either by invoking the distribution of E c , or by distribution 
of a, but the underlying assumption behind the second proposition is hard to justify. We 
have found that the Gaussian distribution of E c i.e. 4>{E C ) = exp[— (E c — E c q) 2 /Aa 2 ] and 
1^(0, T) = /i J °° (f)(E c ) exp(—E c / k,BT)dE c explain the data excellently(shown in Fig.3). E c0 , 
which is the maximum activation energy and a c , are the center and width of the Gaussian 
distribution and found to be 540meV and 80meV, respectively. In this case, the tempera- 
ture dependence of mobility follows ln(fi) oc 1/ksT instead of ln(fi) oc 1/(/c j bT) 2 according 
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to CGDM0|. The spatial correlation in CGDM 7] has been shown to be resulted from 
the long-range interaction between charge carriers and permanent dipole moments of doped 
molecules in polymers and/or host molecules. However, it has been pointed out 9] that 
the mechanism responsible for PF behavior in different conjugated polymers and molecules 
cannot be due to charge-dipole interaction, because the PF behavior of /i has been univer- 
sally observed in several doped and undoped organic semiconductors based on polymers or 
molecules with or without dipole moment. Recently, Yu et. al.(9( have shown using first 
principle quantum chemical calculation that the thermal fluctuations in the molecular ge- 
ometry can lead to spatial correlation through intermolecular restoring force and different 
temperature dependence in this case arises due to the fact that energetic disorder is temper- 
ature independent in CGDM, whereas it increases with temperature in molecular geometry 
fluctuation modei0,|21| resulting ln(fi) oc 1/ksT temperature dependence. We have used 
highly symmetric molecules with negligible dipole moment for our investigations and we 
find it hard to justify strong long-range charge-dipole interaction to be the origin of spatial 
correlation. 

Fig.4 shows the absorption spectra of CuPc. The characteristic two-humped spectra are 
known as Q-band(14) in CuPc. The absorption co-efficient a(E) for a given photon energy 
E is proportional to the probability P{f for the transition from initial state i(in HOMO) 
to the final state /(in LUMO) and the density of the initial and final state and given by, 
a(E) oc T,PifPiPf, where Pi and Pf are the density of states in HOMO and LUMO of CuPc, 
respectively. There are two important parameters involved in the absorption spectra. First 
is the center of the peak at around 2eV, which is the energetic separation between maxima of 
HOMO and LUMO and second is the full width at half maximum of about 460meV, which 
should be the maximum energetic separation between two hopping sites within the intrinsic 
GDOS, as shown schematically in Fig.4. Hence, the maximum contribution in the activation 
energy from static disorder will be around 460meV. We have also found that the 540meV 
is the maximum activation energy (E c o), which has both contributions from static disorder 
and polaron binding energy. Combining these two electrical and optical measurements, we 
found 80 meV(=540-460meV) is the polaron binding energy. If we include the lower and 
upper limits of activation energy due to static disorder, which are 500meV and 580meV 
respectively, the polaron binding energy varies between 40 to 120meV. These observations 
are consistent with the third viewpoint i.e. disordered polaron moving in a correlated energy 



5 



landscape is responsible for charge transport in disordered organic molecular semiconductors. 
These experimental findings settle an extremely important issue regarding the nature of 
charge carriers and have shown conclusively that polarons moving among static disorders 
are responsible for charge carrier transport in disorder molecular semiconductors. 
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Figure captions 

1 J-V characteristics of lOOnm CuPc based single layer device with ITO cathode contact and 

Al anode contact at different temperatures starting at 300K and then at the interval 
of 30K. Experimental data are shown by the symbols and solid lines represent the 
simulated data. Inset shows the log of conductivity of CuPc vs. 1/T 1 / 4 . Empty circles 
are experimental data and solid line is fit with straight line. 

2 Temperature dependence of zero field mobility against 1/T(a) and 1/T 2 (b). Solid lines 

are linear fit to data, which results in (a) /i = 2 x 20~ 4 cm 2 /Vsec, A = A50meV and 
in (b) fi Q = 2 x 10~ 8 cm 2 /V sec, a = 85meV. 

3 Theoretical fit to temperature dependence of zero field mobility in CuPc for two different 

thicknesses using Gaussian distribution of activation energy, which is shown in inset. 

4 Absorption spectra of Cu-Pc for two different thicknesses at room temperature. The 

spectral characteristics and width of the spectra do not depend on the thickness of 
the sample. Inset shows the hopping sites and hopping processes inside the intrinsic 
density of states of disordered molecular semiconductors. 
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Figure.1 
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Figure.3 
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